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Three mode parametric instability has been predicted in Advanced gravitational wave detectors. 

Here we present the first observation of this phenomenon in a large scale suspended optical cavity 
designed to be comparable to those of advanced gravitational wave detectors. Our results show that 
previous modelling assumptions that transverse optical modes are stable in frequency except for 
frequency drifts on a thermal deformation time scale is unlikely to be valid for suspended mass optical 
cavities. We demonstrate that mirror figure errors cause a dependence of transverse mode offset 
frequency on spot position. Combined with low frequency residual motion of suspended mirrors, this 
leads to transverse mode frequency modulation which suppresses the effective parametric gain. We 
show that this gain suppression mechanism can be enhanced by laser spot dithering or fast thermal 
modulation. Using Advanced LIGO test mass data and thermal modelling we show that gain 
suppression factors of 10-20 could be achieved for individual modes, sufficient to greatly ameliorate 
the parametric instability problem. 


I. INTRODUCTION 

Advanced laser interferometer gravitational wave de¬ 
tectors are currently being commissioned ELS- Once 
they reach target sensitivity they have a high probabil¬ 
ity of observing gravitational waves, especially from the 
coalescence of binary neutron stars. Target sensitivity 
requires very high optical power in the detector optical 
cavities, which can allow radiation pressure induced in¬ 
stabilities. 

In 2001, Braginsky et al. SI predicted that opto- 
acoustic interactions in such detectors could lead to a new 
form of instability called a three-mode parametric insta¬ 
bility. It could arise from optical transitions between cav¬ 
ity modes mediated by test mass acoustic modes. Specif¬ 
ically photons from the main interferometer pump mode 
are scattered from thermally excited acoustic modes in 
the test masses. The pump photon creates a phonon- 
photon pair. If the phonon is resonant in a test mass 
acoustic mode, and the photon is resonant in an inter¬ 
ferometer cavity transverse mode, this scattering process 
will occur resonantly. Assuming that the scattered pho¬ 
tons have lower frequency than the pump photons, en¬ 
ergy conservation requires the phonon to increase the oc¬ 
cupation number of the acoustic mode. If the acoustic 
energy injection by this mechanism exceeds the charac¬ 
teristic losses of the acoustic mode, the scattering will 
lead to an exponential growth of the acoustic mode oc¬ 
cupation number. Braginsky showed that the amplitude 
of this scattering process could be large if the spatial 
2D surface amplitude distribution of the acoustic mode 
overlapped the spatial intensity distribution of the optical 
mode, thereby causing three mode parametric instability. 


Subsequently Zhao et al. j5j demonstrated that for 
realistic interferometer designs there was a substantial 
risk of instability, because the high acoustic mode den¬ 
sity in the 50-150kHz range led to numerous accidental 
overlaps of both mode shape and frequency. Such insta¬ 
bilities could not be completely avoided through optical 
design. This led to research focussed on observation and 
study of three mode interactions Q, and on methods for 
suppressing instability [ffl- 

Strigin et a/. El extended the theory to a dual recy¬ 
cling interferometer detector and showed that the multi- 
cavity coupling could reduce the effective linewidth to a 
sub-Hz range. If the high order cavity mode involved in 
parametric instability is resonant in both the arm cavities 
and the recycling cavity, extremely high three mode para¬ 
metric gain could occur. Detailed analysis of a dual recy- 
cling interferometer with realistic test masses by Gras et 
ql.[13] showed that the highest gain could reach ~ 1000 
corresponding to acoustic ring-up times ~ seconds. 

Recently modelling that takes into account large acous¬ 
tic amplitudes, and using parameters close to those of 
Advanced LIGO (aLIGO) has shown that the growth 
of instability saturates. Danilishin et al. showed that 
parametric instability is likely to grow on a time scale of 
minutes for realistic parameters [l4|.fl6j. 

Three mode parametric interactions are extremely sen¬ 
sitive to test mass mirror parameters. This extreme 
sensitivity was emphasised by Ju et al. El who showed 
that mirror radius of curvature changes corresponding to 
wavefront deformations of 10 -6 A could easily be observed 
by monitoring three mode interactions in Advanced in¬ 
terferometers. 

To date three-mode instability has been reported in 
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one free space cavity ex per iment using a picogram mem¬ 
brane in a 10cm cavity [Ij], and in aLIGO jl5l |. 

At the Gingin High Optical Power Facility pj|| a 74m 
optical cavity has been set up to be comparable to the 
conditions of Advanced interferometers. This paper is 
based on observations in this facility which, while demon¬ 
strating instability, has revealed a phenomenon that sup¬ 
presses the exponential growth of instability at low am¬ 
plitudes. 

Previous modelling has ignored two real world aspects 
of practical suspended mass interferometers: a) that mir¬ 
rors after coating have figure errors ~ lnm RMS over the 
central diameter of 160mm, and b) that the laser spot 
position on the mirrors fluctuates due to residual low 
frequency seismic motion. The presence of figure errors 
means that the average radius of curvature of the region 
of the mirror intercepted by the laser beam depends on 
the beam location. This radius of curvature determines 
the transverse mode offset frequency. Because low fre¬ 
quency fluctuations of the spot position causes the laser 
spot to intercept different regions of the mirror surface, 
it follows that there will be dynamical modulation of the 
optical transverse mode offset frequency. The frequency 
modulation causes the parametric gain to be time depen¬ 
dent, and if the modulation amplitude exceeds the trans¬ 
verse mode optical linewidth, the gain can be strongly 
modulated. This can create a situation where paramet¬ 
ric instability does not have time to develop because it 
is only on-resonance intermittently, and for too short a 
time for instability to grow to problem levels. 

In this paper we will show that the above phenomenon 
is likely to reduce the average parametric gain of the can¬ 
didate modes most likely to become unstable, thereby 
significantly reducing the risk of instability. Results are 
confirmed by modelling and by measurements on a 74m 
optical cavity at Gingin. Recognition of this frequency 
modulation suppression mechanism also leads to methods 
by which suppression can be enhanced either by modu¬ 
lated thermal actuation or spot position dithering at fre¬ 
quencies below the gravitational wave sensitivity band. 

In section El we summarise the theory of parametric 
instability and present modelling results showing how 
individual unstable modes can be suppressed by seis¬ 
mic induced frequency modulation. In section IIIII we 
use aLIGO test mass mirror metrology data to estimate 
the frequency modulation expected for small spot po¬ 
sition motions in aLIGO. In section IIVI we present re¬ 
sults obtained with the Gingin high optical power cav¬ 
ity: both the observation of parametric instability and 
the frequency modulation that greatly reduces the risk 
of instability. We discuss the results obtained, and their 
implications for aLIGO. We also present thermal actua¬ 
tion modelling results to estimate the suppression factors 
achievable. 


II. THEORY OF PI AND EFFECT OF 
TRANSVERSE MODE FREQUENCY 
MODULATION 

Three mode opto-acoustic interactions occur when the 
frequency difference between an optical cavity pump 
mode at frequency wo and a transverse mode at frequency 
wi is appropriately tuned to the frequency of an acoustic 
mode at frequency A m . This three mode interaction reso¬ 
nance is defined by A m = (wo — wi) — w m = Aw — w m = 0. 
The parametric gain R characterises the ratio of acoustic 
energy input compared to mirror acoustic mode losses. 
If R > 1, the system is acoustically unstable, and the 
acoustic mode will grow exponentially until either non- 
linearities cause saturation [h|, or else the cavity loses 
lock. In this paper we are concerned only with small am¬ 
plitude excitation so can ignore non-linearities. The mag¬ 
nitude of R depends on cavity input power, on acoustic 
and optical mode losses, and on the spatial overlap be¬ 
tween the relevant modes. For any pair of acoustic and 
optical modes, the gain R can be expressed as[3] 

PAloi 1 
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Here P is the input power to the cavity, 70 , 71 , and 
7 m are the half-linewidth of the two optical modes and 
the acoustic mode of the test mass respectively, M is the 
mass of the test mass, L is the length of the cavity. The 
A is the overlap factor including the mass to the effective 
mass ratio as defined in Q. The optical mode spacing 
Aw is a function of the radius of curvature of the mirrors 
of the optical cavity and is given by 

Aw = |-(m + n) cos- 1 ^- - -7-)), (2) 

where i?i and R 2 are the radii of curvature of the end 
mirrors of the cavity, m and n are integers describing 
the order of the optical mode. The sign depends on 
the cavity configuration. Equation @ assumes perfect 
spherical mirrors, but we will assume that in the case of 
figure errors the mode spacing is defined by the average 
radius of curvature at the laser spot position, averaged 
over the effective spot size. 

Equation m considers only the Stokes process where 
parametric amplification or instability processes occur 
due to a single high order optical mode. Here we want to 
focus particularly on the case where dynamic detuning 
causes A m to be time dependent. We consider the case 
of harmonic detuning given by 

A m (f) = AmOCOSW^t, (3) 

where w^ is a dynamic tuning frequency. In suspended 
mass interferometers the test mass-mirrors are supported 
by low frequency pendula which isolate against vibration. 
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The test mass positions are controlled by feedback, but 
finite residual motion is inevitable because of the require¬ 
ment that the test masses be inertial within the gravita¬ 
tional wave signal band. 

Thus in practice, test masses can be expected to have 
significant motion at pendulum normal mode frequencies 
0.1-lHz. This gives rise to a modulation in the laser spot 
position. If the mirrors are imperfect, the mirror radius 
of curvature (averaged over the laser spot size) will vary 
smoothly with spot position. In this case, modulation 
in spot position can modulate the transverse mode offset 
frequency, thus causing time dependent detuning fluctu¬ 
ations. 

Spot position motion will also modulate the modal 
overlap parameter. However for millimetre scale mo¬ 
tions the overlap parameter modulation is small com¬ 
pared with the effect of detuning, and is ignored in the 
following analysis. 

Assuming that Aw m changes according to Eq ([3]), the 
parametric gain is given by 


m 


R 


max 


1 + (a cosUdt) 2 ’ 


(4) 


where a = A m o/7i is the normalised frequency detun¬ 
ing modulation amplitude. Equation [4] allows estimation 
of the effects of modulation on the growth of paramet¬ 
ric instability. As discussed above, modelling has shown 
that the characteristic ring up time scale for parametric 
instability in a detector similar to aLIGO is likely to be 
~ 10 2 s jjj|. Since u>d is fast compared with such ring 
up times, one would expect to observe modulated signal 
growth. 
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FIG. 2. Effective parametric gain suppression as a function 
of dynamic detuning amplitude a. 


a typical acoustic mode decay time r without need to 
specify acoustic mode frequency. We chose to use r = 6s, 
corresponding to acoustic quality factor Q m = 10 6 and 
4 x 10 6 for frequencies 50kHz and 200kHz respectively. 
Results are compared with acoustic mode ring up curves 
in the absence of dynamic detuning (a = 0) for R = 6. 
In the case of a = 4 the ring up slope is equivalent to 
that of a system with a = 0 and R = 1.45 as indicated in 
the figure. This represents a gain suppression factor ~ 4. 
Clearly in all cases, frequency modulation suppresses the 
effective parametric gain as determined by the average 
slope of the ring up curves. For a = 10 we see that insta¬ 
bility has been replaced by a modulated acoustic mode 
amplitude which while not harmonic, is stable in time. 
The equivalent parametric gain R a in the presence of 
harmonic dynamic detuning is given by: 



FIG. 1. Acoustic mode amplitude ring up curves for var¬ 
ious detuning amplitudes. Here we assume maximum gain 
Rmax= 6 , an acoustic mode ring down time r= 6 and a dy¬ 
namic modulation frequency of 0.1 Hz. For comparison, cases 
for on resonance (a=0) with R ma x= 6 and R m ax= 1-45 are also 
plotted. 

Figure [T| shows examples of possible acoustic mode 
ring up signatures. We assume parameters compara¬ 
ble to those of the experiment reported in this paper: 
fd = w<i/27 t = 0.1Hz, Rmax = 6, and a normalised detun¬ 
ing amplitude a = 1, 2, 4 and 10. It is sufficient to choose 


R 
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The suppression of effective parametric gain as a func¬ 
tion of modulation amplitude is shown in Fig. [2] Para¬ 
metric gain can be suppressed by an order of magnitude 
for a = 10. 

The mechanism discussed above occurs because the dy¬ 
namic detuning modulation frequency is fast compared 
with the acoustic mode ring up time scale. The observed 
acoustic mode amplitude modulation occurs at double 
the dynamic detuning frequency u>d- While the effective 
parametric gain is independent of u>d, the peak to peak 
acoustic mode amplitude within one cycle is inversely de¬ 
pendent on u>d- Figure [3] shows some examples for three 
different dynamic detuning frequencies. The amplitude 
modulation waveform is highly non-linear since it is due 
to a Lorentzian modulation acting on the exponent. 

It can be seen that if the detuning frequency is too 
slow, the acoustic mode amplitude can grow to a very 
large value within half a detuning period. It is possible 
to define a lower limit for the dynamic detuning frequency 
fdjim to prevent the acoustic mode amplitude excursion 
from exceeding /? times its original value within one cy¬ 
cle. Figure 0] shows three curves showing the lower limit 
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FIG. 3. As the dynamic detuning frequency increases the 
acoustic mode amplitude excursions are reduced but the ef¬ 
fective parametric gain is unaltered. Here three detuning fre¬ 
quencies 1Hz, 0.1Hz and 0.01Hz are shown. The detuning 
amplitude is fixed at a = 5 (Rmax = 6 and r = 6). 


of the dynamic detuning frequency as a function of detun¬ 
ing amplitude a, for two values of R m ax and two values 
of the acoustic amplitude excursion limit /3. For example 
if Rmax — 10 , and r = 6 s, and a requirement /3 = 2 , 
then the dynamic detuning frequency is limited within 
the range 0.1Hz - 0.6Hz assuming detuning modulation 
amplitudes a is between 2 and 16. We see that larger 
detuning frequencies or larger detuning amplitudes both 
act to reduce amplitude excursions. This defines the pa¬ 
rameter space for suppressing parametric instability by 
the dynamic detuning mechanism. In the next sections 
we will see that the above mechanism can occur naturally 
as a result of residual motion in the presence of mirror 
figure errors, which we consider in the context of aLIGO 
and a 74m cavity at Gingin. 
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FIG. 4. Dynamic detuning frequency limit to prevent am¬ 
plitude excursions exceeding a predetermined value /?■ For 
example, for Rmax = 10 and r = 6s, with an amplitude 
growth requirement of /3 = 2, then the minimum dynamic 
detuning frequencies is limited to between 0.6 Hz and 0.1 Hz 
for detuning amplitude a in the range of 2 and 16. 


III. FREQUENCY MODULATION BY MIRROR 
FIGURE ERRORS 

The test mass mirrors in aLIGO have radii of curva¬ 
ture of ~ 2000m. Figure errors mean that the effective 
radius of curvature depends on the spot position. For 
example, using equation [2j with aLIGO arm cavities, if 
the ETM RoC changes lm from nominal value of 2242m 
(corresponding to a sagitta change within a beam diam¬ 
eter ~ 0.3nm), the TEM10 cavity mode frequency will 
change ~ 13Hz. 

The residual motion of the test masses in interferom¬ 
eter arm cavities therefore cause high order mode fre¬ 
quency modulation. Residual angular motion creates 
beam residual motion on the test mass surface mil¬ 
limetres 0. Depending on the test mass figure errors, 
this residual motion causes dynamic detuning of the cav¬ 
ity high order mode frequency at the frequencies of test 
mass pitch and yaw motion. 



FIG. 5. a) Typical aLIGO test mass figure errors (compared 
to a perfect sphere of RoC ~ 2242m) showing deformations 
across a mirror diameter, b) FFT code model for frequency 
offset as a function of test mass angular motion. 

To estimate the aLIGO arm cavity high order mode 
frequency changes as a function of the residual test mass 
ang ular motion, we used measured test mass surface data 
[20| in an interferometer simulation codes (OSCAR [22j 
and FOPGUI) to simulate cavity transverse mode de¬ 
tuning. For the simulation we fixed the ITM and mod¬ 
elled ETM misalignment at various angles from zero to 
0.25 micro-radians. Figure[5Ka) shows an example of the 
input data in the form of a cross section across a test mass 
diameter showing how the figure errors increase with ra¬ 
dius. Figure [5Kb) shows the calculated data for real 2D 
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surface profiles. 

Figure [5] (b) shows that the cavity mode spacing in¬ 
creases roughly quadratically with ETM misalignment 
angle. Note that 0.1 microradian corresponding to 
~2mm beam position displacement on the test mass in 
a typical aLIGO arm cavity. We extended the simula¬ 
tion to the test mass rotation corresponding to a spot 
displacement ^6 mm. We note that in reality both test 
masses move independently of each other, thereby creat¬ 
ing somewhat larger detuning amplitudes. 

Cavity high order mode frequency modulation can also 
be artificially created by applying modulated heating to 
the test mass. We used ANSYS FEM simulation software 
package to simulate the transient thermal deformation of 
the test mass surface under sinusoid heating power. Fig¬ 
ure [6] shows the maximum thermal deformation when a 
0.1 Hz modulated heating beam of 50mm radius with 2W 
peak-to-peak power amplitude is applied to the test mass 
front surface. This deformation corresponds to a cavity 
mode spacing frequency change of ~40Hz in aLIGO arm 
cavity simulated using FFT code [23]. 



y (mm) - 200 -200 x(mm ) 


FIG. 6. The maximum thermal deformation when 0.1 Hz 
sinusoid heating power of amplitude 2W applied on the front 
surface of the test mass 

The above results indicate that passive detuning fre¬ 
quency modulation in aLIGO would be expected to be 
~few Hz for the TEM10 transverse mode, which is much 
smaller than the arm cavity linewidth and does not have 
significant effect on parametric gain. However this could 
be increased to ~40Hz modulation using C02 laser heat¬ 
ing. It is important to note that the highest predicted 
parametric gains in aLIGO are for modes up to 4th or¬ 
der Q. By equation [2] detuning scales with mode order. 
Thus the above estimates would correspond to at least 
4 times larger modulation (~160Hz) for 4th order insta¬ 
bilities. The aLIGO arm cavity half-linewidth is ~40Hz. 
Thus the parametric gain associated with arm cavity op¬ 
tical modes would be suppressed by factors ~few, and for 
low order cavity modes the above modulations could be 
negligible. Reference @ shows that the highest paramet¬ 
ric gain instabilities are associated with modes that are 
resonant in the power recycling cavity. For these modes, 
the coupled cavity linewidth is about 0.3 Hz [u|, the 
normalised dynamic detuning amplitude can exceed 50 


times of the coupled cavity linewidth. Thus it is most 
likely that intrinsic passive detuning in aLIGO will lead 
to a parametric gain suppression factor > 100 for those 
modes resonant inside the recycling cavity. A detailed 
simulation to explore how the cavity high order mode 
frequency modulation affects the broad spectrum of para¬ 
metric instability in aLIGO is beyond the scope of this 
paper. However the experimental observation of para¬ 
metric instability presented in the next section largely 
confirms the above theory. 


IV. HIGH OPTICAL POWER CAVITY 
OBSERVATIONS 

We studied three mode parametric instabilities at the 
Gingin high optical power facility. The experimental 
setup is shown in Figure [7] A 74m long optical cavity 
with fused silica test masses is suspended from high per¬ 
formance vibration isolators [23l . |24I | by a modular 4-wire 
test mass suspension system developed at UWA. The test 
masses are installed in two large vacuum chambers con¬ 
nected by a 400mm vacuum pipe. The system was as¬ 
sembled in clean room conditions and uses a hydrocarbon 
free vacuum system to enable high optical power densities 
to be achieved. 

Both test masses are 50nnn in diameter and 50mm 
thick, with mass ~0.8kg. The nominal RoC of the two 
test masses are 37.5m and 37.4m. The test masses have 
a very sparse mode spectrum compared to aLIGO test 
masses, so that three mode parametric interactions need 
to be tuned to specific candidate acoustic modes. This is 
achieved by using a power stabilised C02 laser to ther¬ 
mally tune the ITM RoC to create three mode tuning for 
the specific candidate acoustic modes [25| . The dominant 
test mass residual angular motions are at frequencies of 
0.15 Hz. 

The measured cavity finess is 14500±300. The light 
source is a 50W fibre laser amplifier fed by a 400mW 
Nd:YAG NPRO seed laser. The seed laser is frequency 
locked to the long cavity using PDH locking [26]. The 
cavity transmission is detected by a quadrant photodiode 
(QPD). The differential output of the QPD measures the 
beating between the cavity fundamental mode and the 
first order mode while the sum of the QPD output mea¬ 
sures the total cavity transmitted power. A spectrum 
analyser (Agilent 89410A) and a PC are used to analyse 
and to record the signal. 

Using the ANSYS software package, we first analysed 
the test mass acoustic mode structure and frequencies. 
Based on the simulation we then identified one particu¬ 
lar acoustic mode that has good overlap with the cavity 
first order mode, and minimum vibration amplitude at 
the suspension point to minimise the mechanical loss in¬ 
troduced by the suspension. Our target mode, with sim¬ 
ulation frequency 150.49 kHz is in the range for easy CO 2 
laser thermal tuning. The mode amplitude distribution 
on the test mass surface is shown in Fig. [8] The over- 
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Fibre Laser 



FIG. 7. Schematic diagram of the experimental setup: The 
laser light from a seed laser is amplified by a 50 W fibre laser 
amplifier. The high optical power laser beam is injected into 
the 74m long optical cavity. The seed laser is frequency locked 
to the long cavity using PDH locking. The cavity transmitted 
beam is detected by a quadrant photodiode (QPD). The dif¬ 
ferential signal from the QPD measures the beating between 
the cavity fundamental mode and the first order mode. 



FIG. 8. The test mode acoustic mode amplitude distribution 
on the surface. 


lap factor taking into account the mode effective mass is 
~16. The measured mode frequency is ^150.2 kHz (de¬ 
pending on the temperature.) The measured mechanical 
Q-factor using the ring-down method is ~ 3.4 x 10 6 . 

Three mode interaction conditions are achieved by 
tuning the TEM00 and TEM10 mode spacing close to 
150.2kHz using C02 laser thermal tuning. Measure¬ 
ment of the tuning is relatively easy because residual 
laser beam-jitter noise gives rise to a small amount of 
TEM10 mode power inside the cavity which beats with 
the TEM00 at the QPD, allowing the TEM10 offset fre¬ 
quency to be monitored as a beat note. This provides 
a means for monitoring the mode spacing by measuring 
the cavity transmitted power on the QPD where the two 
modes are mixed. 

The mode spacing was observed to fluctuate with a 
typical peak to peak amplitude ~ few kHz. To con¬ 
firm that these fluctuations were associated with beam 
spot position on the test masses, we recorded the cavity 
mode spacing and the beam position simultaneously for 
the ITM. Figure [9] shows the mode spacing as a function 
of the beam position on the ITM in horizontal direc¬ 
tion. The beam position was determined by recording 
the video of the CCD camera and then analysed by ref¬ 



Horizontal spot position (nun) 

FIG. 9. The correlation between horizontal spot position of 
the ITM alone and the transverse mode frequency indicates 
that the frequency detuning is caused by the spot position 
change. The solid line is the linear fit to the measurement 
data. The spread of the data is due to the fact that the laser 
spot position is also change in ETM. 


erencing it to the test mass diameter. In the horizontal 
direction, there is a linear correlation between increased 
mode spacing with increased beam position. The solid 
line in figure [5] is a linear least squares fitting to the mea¬ 
surement data. The relative large scatter is due to the 
fact that we recorded only the ITM beam position while 
the ETM beam position is also not stable. The effect is 
more difficult to measure in the other axis because the 
suspensions introduce much smaller vertical beam posi¬ 
tion fluctuations. However the single axis correlation is 
sufficient to confirm our conjecture that mirror figure er¬ 
rors translate into dynamic detuning. 

We do not have precise metrology of our test mass mir¬ 
ror profiles. However the observed fluctuations are con¬ 
sistent with the mirror figure error specification of lnm. 
It is interesting to note that in principle simultaneous 
measurement of spot position on both test masses and 
transverse mode frequency offset could be used to allow 
precise metrology of both test masses. 

When the cavity is correctly tuned the three mode in¬ 
teraction occurs, and the signal at the QPD becomes 
dominated by the beating between TEM00 and TEM10 
modes at the acoustic mode frequency. The signal is pro¬ 
portional to the acoustic mode amplitude, the TEM00 
mode power and the TEM10 mode detuning. The signal 
is normally most easily observed by mixing the acous¬ 
tic frequency with a local oscillator, combined with a low 
pass filter, so as to reduce the signal frequency to < 10Hz. 

As discussed above residual motion causes cavity de¬ 
tuning. The residual motion amplitude depends on en¬ 
vironmental noise, which excites the suspension normal 
modes. Most of the time we observe dynamic detun¬ 
ing with a frequency amplitude of 1-5 kHz. Even un¬ 
der these circumstances the acoustic mode signal at fre¬ 
quency ~150.2 kHz can normally be clearly observed. 

Wind forces on the laboratory, microseismic activity 
and human activity all contribute to degrading the resid¬ 
ual motion. During quiet times the residual motion is 
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FIG. 10. The QPD differential output signal at test mass 
acoustic mode frequency (150.28 kHz). The signal was down- 
converted to ~0.91 Hz by mixing with an local oscillator sig¬ 
nal. The solid line is a fitting curve of 0.91Hz with para¬ 
metric gain Rmax = 6 and detuning amplitude a=2. The 
growing signal envelope (dashed line) is consistent with sus¬ 
pension modulation at 0.15Hz. The effective parametric gain 
is ~1.45. 


reduced and for times of ~30 seconds the detuning am¬ 
plitudes can be less than a few cavity linewidths. In these 
short periods of time conditions are suitable for observing 
three mode parametric instability. 

To observe the signature of parametric instability we 
increased the cavity circulating power to ~ 30 kW. For 
periods of time ~10 to 30 seconds, when the dynamic 
detuning is low, the acoustic signal can be observed ring¬ 
ing up with time, as shown in figure flOl In this case the 
acoustic signal frequency was down converted to 0.91 Hz 
as discussed above. 

Observations under best tuned quiet conditions show 
the acoustic signal growing for times ~ 14 seconds. The 
amplitude growth is modulated but more complex than 
the single modulation frequency model used in section 
IE due to the presence of several low frequency modula¬ 
tion frequencies associated with the angular motion both 
test masses. Beating also occurs, due to the fact that 
the two test masses have closely spaced suspension nor¬ 
mal modes. The beating causes the detuning amplitude 
to vary periodically over time scales ~ 30 seconds. The 
effective parametric gain based on observed ring-ups dur¬ 
ing times of minimum detuning amplitude such as shown 


in Figure [10] is R ~1.45. 

In Figure [El we have fitted to a double frequency of 
a single 0.15Hz suspension mode to model the dynamic 
detuning. This gives a modest fit to the data but a com¬ 
plete fit is not possible due to the stochastic nature of 
the seismic excitation of the normal modes. 


V. CONCLUSIONS 

We have created conditions in which three mode para¬ 
metric instability can occur in a suspended high power 
optical cavity designed to mimic conditions comparable 
to those expected in advanced gravitational wave detec¬ 
tors. We have observed time dependent growth of a 
150.2kHz acoustic mode, consistent with a new model 
of parametric instability for suspended mass optical cav¬ 
ities. The gain in the parametric instability regime is 
lower than previously expected, and modulated by low 
frequency residual motion. Results are consistent with 
a new model for the build up of instability in which 
transverse mode frequency fluctuations act to reduce the 
parametric instability power build up through dynamic 
detuning which itself is caused by residual motion in the 
presence of nm-level mirror figure errors. Data on aLIGO 
optical cavities indicate that the same phenomenon will 
act to reduce the risk of parametric instability for the 
highest parametric gain modes. Mirror imperfections 
have beneficial effects in this regard. Results also point to 
simple methods for reducing parametric gain by thermal 
modulation of test masses or by low frequency dithering 
of the test masses. Further studies on full scale detec¬ 
tors to quantify the dynamic detuning and linewidths of 
transverse modes are needed to quantify these effects. 
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